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Neutron scattering measurements have been performed on polycrystalline samples of the newly 
discovered layered superconductor LaOo.5Fo.5BiS2, and its nonsuperconducting parent compound 
LaOBiS2- The crystal structures and vibrational modes have been examined. Upon F-doping, 
while the lattice contracts significantly along c and expands slightly along a, the buckling of the 
BiS2 plane remains almost the same. In the inelastic measurements, a large difference in the high 
energy phonon modes was observed upon F substitution. Alternatively, the low energy modes 
remain almost unchanged between non-superconducting and superconducting states either by F- 
doping or by cooling through the transition temperature. Using density functional perturbation 
theory we identify the phonon modes, and estimate the phonon density of states. We compare these 
calculations to the current measurements and other theoretical studies of this new superconducting 
material. 

PACS numbers: 61.05.F-,63.20.dk,74.25.Kc,78.70.Nx 



Superconductivity has been fascinating scientists for 
more than a century [1]. Several families of supercon- 
ducting (SC) materials and mechanisms have been found 
and proposed: notably type I or BCS superconductors 
where superconductivity is mediated by phonons, i.e. lat- 
tice vibrations, and type II superconductors whose su- 
perconductivity is yet to be fully understood. Many ma- 
terials that belong to the latter category have layered 
crystal structures with low dimensionality. The cuprates 
with Cu02 layers and the Fe-based superconductors with 
Fe-An (An: pnictogen or chalcogen anion) layers are 
two extensively studied examples 0-0] • Very recently, 
a new family of materials based on BiS2 layers has been 
found to be SC at low temperatures: Bi^^SO^i-^ [H, 
LaOi_,F,BiS 2 @, NdOi^F.BiSs 0, PrOi_,F,BiS 2 
|H, and CeOi_^F cc BiS2 Q. The natural question that 
arises is whether the new Bi-based superconductors are 
type I BCS superconductors or another family of type II 
superconductors yielding a new route for unconventional 
superconductivity. 

Several theoretical studies have been reported espe- 
cially for LaOo.5F .5BiS 2 with the highest T c w 10.8 
K among the Bi-superconductors until now. Theoreti- 
cal consensus so far is that in LaOo.5Fo.5BiS2 the Fermi 
level crosses conduction bands yielding electron pock- 
ets [TMl- 

From the quasi-one-dimensional (Id) nature 
of the conduction bands, these Fermi surfaces(FS) are 
nested by a (7r,7r) wave vector. The many body inter- 
actions and their relation to the superconductivity are, 
however, still controversial. A group of theories propose 
that FS nesting, just as in the Fe-pnictides, makes elec- 



tronic correlations a candidate for the pairing mechanism 
111 . [Til . Alternatively, electron-phonon (e-ph) coupling 
constant (A) calculations 13, 3, lfij find A ~ 0.85 with 
T c « 11.3 K, close to the experimentally found value. 
This suggests LaOi_^F :E BiS2 can be a conventional su- 
perconductor with strong e-ph coup ling . Moreover, the- 
ories predict structural instabilities 13, which places 
this system in close proximity to competing ferroelectric 
and charge density wave (CDW) phases [lfj. This insta- 
bility is believed to be due to the anharmonic potential of 
the S ions lying on the same plane as the Bi ions. Just like 
the spin-density wave (SDW) instability in Fe-based su- 
perconductors, a CDW instability from negative phonon 
modes at or around the M point, (7r,7r), is suggested to be 
essential to the superconductivity. The buckling of the 
BiS2 plane, being related to the electronic structure and 
FS, was also suggested to be crucial in inducing supercon- 
ductivity, and expected to decrease in the SC compound 

Experimentally, doping dependent structural measure- 
ments [17j suggest that the a lattice parameter is relevant 
to superconductivity. Both a and T c increase with F 
substitution, and become maximal near optimal doping 
x w 0.5. The existence of electron pockets near the FS 
was experimentally shown by Hall effect measurements 
18l LL9j. Other experimental studies to test the theo- 
retical scenarios, however, are limited. Thus, detailed 
experimental studies of the crystal structure and lattice 
vibrations of non-SC and SC compounds are crucial in 
understanding the new Bi-superconductors. 

Here we report neutron diffraction and inelastic neu- 
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FIG. 1. (color online). Temperature dependent resistivity of 
LaOBiS2 (black lines) and LaOo.5Fo.5BiS2 (red squares) . The 
inset shows the crystal structure of LaOi- cc F cc BiS2. 



tron scattering measurements both on non-SC LaOBiS2 
and SC LaOo.5Fo.5BiS2 to investigate how the crystal 
structure and phonon density of states (PDOS) change 
as a function of F-doping and temperature. We also per- 
formed density functional linear response calculations to 
identify the corresponding phonon modes. Experimen- 
tally obtained structural parameters are compared be- 
tween the non-SC and SC compound and with calculated 
values. While significant differences in the phonon spec- 
trum was observed at higher energy transfers upon F- 
substitution, little change is observed in the low energy 
portion of the spectrum. The calculations suggest sig- 
nificant modification in the low energy spectrum which 
should be relevant to the superconductivity, contrary 
to the experiment. We do not observe any meaningful 
change as a function of temperature in the phonon spec- 
trum of the SC compound in the vicinity of T c . 

A 4.89 g polycrystalline sample of LaOBiS2 was syn- 
thesized using the solid state reaction method under am- 
bient pressure while a 0.89 g polycrystalline sample of 
LaOo.5Fo.5BiS2 was synthesized under high pressure at 
the NIMS, Tsukuba, Japan. Details concerning sample 
synthesis are described in Ref. 0, [ItJ . Figure [T] shows 
the low-temperature resistivity of the parent and F sub- 
stituted compound. LaOo.5Fo.5BiS2 has a clear SC tran- 
sition at T c « 10.8 K. 

Neutron scattering measurements were performed at 
the Spallation Neutron Source (SNS) using the POW- 
GEN diffractometer, the wide angular range chopper 
spectrometer (ARCS), and the Cold Neutron Chop- 
per Spectrometer (CNCS). Samples were loaded into 
vanadium (at POWGEN and CNCS), or aluminum (at 
ARCS) cans with a He atmosphere and mounted to the 
closed cycle refrigerator (at POWGEN and ARCS) or 
liquid helium cryostat (at CNCS). Neutron diffraction 
data were collected using a wavelength band to cover a 
wide range of d spacing from 0.55 to 4.12 A at POW- 
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FIG. 2. (color online). Neutron powder diffraction data ob- 
tained from (a) LaOBiS2, and (b) LaOo.5Fo.sBiS2 at 15K. 
Black crosses represent observed data. Red, green, and blue 
solid lines are the calculated intensity, estimated background, 
and difference between the observed and calculated intensity 
obtained by GS AS [3 , respectively. 



GEN 20} . Inelastic neutron scattering (INS) measure- 
ments were performed at ARCS [2l| with monochromatic 
neutrons of incident energies Ei = 40 and 80 meV. For 
improved resolution at low energy transfer, INS measure- 
ments were performed at CNCS [22| with E$ = 25 meV. 
All the INS data presented here are corrected for back- 
ground by subtraction of an empty can measurement. 

The phonon calculations were performed using 
QUANTUM-ESPRESSO [23]. A 9 x 9 x 3 regular grid 
over the irreducible Brillouin Zone was used for the self- 
consistent calculation of the F substituted and parent 
compound. A manual check of convergence for grid den- 
sity, energy cutoff, and lattice parameter values was per- 
formed before optimization of structural and atomic po- 
sition values. To simulate the half doping in the SC sam- 
ple, we replace oxygen at one of the 2a Wyckoff sites with 
F in an ordered fashion. 

Figure [2ja) shows neutron diffraction data and Ri- 
etveld refinements obtained from LaOBiS2- The nuclear 
Bragg peaks are instrumental resolution limited indicat- 
ing good crystallinity. The data for the parent compound 
are well reproduced by the PA/nmm space group with 
structural parameters and goodness of fit listed in Ta- 
ble H We determine the crystal structure, as shown in 
the inset of Fig. [TJ to be very similar to that foun d by 
the preliminary x-ray scattering measurements 0, \u\ . 
The Bi ions form a square lattice just as the Cu and Fe 





Experiment 


Calculation 


LaOBiS2 


wRp = 0.0799 




P4/nmm T=15K 


X 2 = 13.48 




a (A) 


4.05735(5) 


4.03949 


c(A) 


13.8402(3) 


14.30361 


La 2c (0.5, 0, z) 


0.09065(9) 


0.08527 


Bi 2c (0, 0.5, z) 


0.3688(1) 


0.36652 


Si 2c (0.5, 0, z) 


0.3836(3) 


0.39493 


D2 2c (U.o, U, zj 


U.81Ul(zj 
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ZSi — Z51 /zBi 


4.Ul^oJ \/o) 




LaO .5F .5BiS 2 


wRp = 0.0647 




P4/nmm T=15K 


X 2 = 2.252 




a (A) 


4.0651(3) 


4.07989 


c(A) 


13.293(7) 


13.42520 


La 2c (0.5, 0, z) 


0.1007(5) 


0.10340 


Bi 2c (0, 0.5, z) 


0.3793(5) 


0.38386 


bi 2c (0.5, 0, zj 


0.362(2) 


0.38401 


S 2 2c (0.5, 0, z) 


0.815(1) 


0.81600 


2a (0, 0, 0) 






F 2a (0, 0, 0) 






\%Bi — ZSl \ /ZBi 


4.5(5) (%) 


0.02413 (%) 



TABLE I. Refined structural parameters of LaOBiS2 and 
LaOo.5Fo.5BiS2 obtained from neutron powder diffraction 
and calculated with structural optimization implemented in 
QUANTUM ESPRESSO. The quantity \z B i - z S i\ /z B i char- 
acterizes the amount of buckling. wRp and x 2 respectively 
are the weighted R factor and chi squared values from the 
structural refinements [24| . 

ions do in the cuprates and Fe-based superconductors, 
respectively. There are two distinct Wyckoff sites for the 
sulfur ions, Si and S2. While Si ions reside nearly on the 
same plane as the Bi square lattice with some buckling, 
S2 ions are just above or below the Bi ions. 

On the other hand, for the SC LaOo.sFo.5BiS2 most of 
the nuclear Bragg peaks are broader than instrumental 
resolution (Fig. Efb)), indicating imperfect crystallinity. 
Similar broadening has been reported in prior X-ray mea- 
surements 0, We note that the broad peaks have 
an asymmetric line-shape, characteristic of a low dimen- 
sional crystal ordering. The asymmetric broad peaks in- 
dex with nonzero I values, while I = Bragg peaks are 
considerably sharper and more symmetric. This is clearly 
seen in Fig. [2fb), for example, where the (200) peak at 
d = 2.03 A and the (110) peak at d = 2.87 A are sharp 
while the (114) peak at d = 2.17 A and the (102) peak 
at d = 3.47 A are broad and asymmetric. This suggests 
that strain may be induced along the c-axis due to a ran- 
dom replacement of the F ions at O sites. As a result, the 
La(0,F) planes are not well ordered along the c-axis. The 
solid line in Fig. [2fb) is our best refinement to the data, 
where a phenomenological model of anisotropic broad- 
ening has been used [25[. This coarsely reproduces the 
diffraction data and allows for reasonable determination 
of the structural parameters, summarized in Table HI 

We find upon substitution with F, the lattice elongates 
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FIG. 3. (color online). GDOS as a function of Huj for the (a) 
parent non-SC and (b) F substituted SC compound. Black 
symbols are the experimental data; blue solid lines are the 
calculated GDOS scaled for comparison. Three measurements 
have been combined: below 10 meV are from CNCS Ei = 25 
meV, between 10 and 30 meV from ARCS E* = 40 meV, 
and above 30 meV from ARCS Ei = 80 meV. To obtain the 
GDOS, the neutron scattering intensity was averaged over 
momentum transfers of [3, 6] A , [4, 7] A , and [3, 8.5] A , 
respectively. The phonon modes at the T point are shown in 
the inset with their symmetry for the parent (above) and SC 
(below) compound. Atoms are shown with the same colors as 
in Fig.[T] and their relative displacements are represented with 
thick blue arrows. The left four figures in the inset are the 
vibrational modes that are theoretically expected to have the 
large e-ph couplings, with each A shown in parentheses. For 
the first mode shown here, the corresponding energy becomes 
negative (—2.9 meV) in the parent compound. For doublet 
modes, only one of the two orthogonal modes is shown. 



along a by « 0.2% while it contracts along c by « 4.1%, 
consistent with previous x-ray results [17| . This F-doping 
induced lattice change is reproduced by our calculations 
although the rates are different. The buckling of the BiS2 
plane is found to slightly increase as shown in Table HI 
This contradicts a prior theoretical prediction [l6| as well 
as our own calculations. 

We note that there is discrepancy between experimen- 
tal and calculated structural parameters. The discrep- 
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FIG. 4. (color online). Temperature dependence of the 
dynamical susceptibility. The data have been measured at 
ARCS with Ei = 40 meV, and integrated over Q = [4, 7] 

A" 1 . 

ancy in the z position of the Si atom was regarded as 
a sign of a possible structural instability in this system 
[l6| . This deviation becomes larger with F substitution. 
According to frozen phonon calculations, the instability 
at or near the M point can reduce the symmetry of the F 
substituted compound statically to P222i [l6[, the CDW 
phase. This structural transition, however, could not be 
experimentally confirmed due to the Bragg peak broad- 
ening. 

To examine if the pairing mechanism of the supercon- 
ductivity is phononic, we have performed INS measure- 
ments to obtain the neutron weighted generalized PDOS 
(GDOS). Figure [3] shows the GDOS as a function of en- 
ergy transfer, Hlj, for the two compounds at T = 5 K. For 
the non-SC LaOBiS2, (Fig. [3fa)), there are well-defined 
phonon modes over a wide range of Hlj up to 70 meV. 
At least five prominent bands of lattice vibrations are 
present at 7.7(2), 14.2(1), 22.9(2), 40.4(1), and 61.8(4) 
meV. Upon F-substitution, as shown in Fig. [3] (b), the 
phonon modes at higher energies change significantly. All 
the vibrational modes become broader than their corre- 
sponding modes of LaOBiS2- Furthermore, the top of the 
band softens in energy from 61.8(4) to 55.2(1) meV, and 
the sharp 40 meV peak significantly weakens and broad- 
ens. Conversely, the first two low energy peaks remain 
similar to those of the parent compound even though 
they broaden. 

Our calculated GDOS are shown as solid lines in Fig. 03 
For the non-SC sample, it reproduces the observed promi- 
nent phonon modes reasonably well. Our calculations at 
the zone center show that the 40 and 62 meV bands are 
mainly due to vibrations of light ions. For example, in 
the insets of Fig. [3j the last 5 figures show the highest 
energy vibration modes. These modes mostly involve O 
and/or S2 vibrations. The intermediate energy modes 
around 13 and 23 meV are mainly due to Si and/or S2 
vibrations, while the low energy modes below « 10 meV 
are due to the vibrations of the BiS2 and/or LaO layers. 



Our calculations at the zone center also yield imaginary 
unstable phonon modes involving vibrations of the BiS2 
layer, which was previously reported to be the sign of 
anharmonic ferroelectric soft phonons [l6| . 

For the SC-sample, the calculated GDOS reproduces 
the features occurring at large huo (albeit shifted in en- 
ergy); the broadening of high energy peaks, reduction 
of the sharp 40 meV peak, and softening of the high- 
est energy O/F vibration modes. These changes can be 
partially understood as being due to the high energy O 
modes being shifted because of substitution with heav- 
ier F ions. This, however, fails to explain the low en- 
ergy data. According to the calculation, most of the 
e-ph coupling comes from the low energy modes below 
« 20 meV, and as such this is where we expect meaning- 
ful changes of the GDOS relevant to the superconduc- 
tivity to occur. While the theory predicts considerable 
re-distribution of GDOS in this low energy region, we do 
not observe any such change in the measured spectrum 
upon F-substitution. 

To further probe for a phonon anomaly associated with 
superconductivity, we also examine the temperature de- 
pendence of the dynamical susceptibility, x ( C<J )? °f the 
SC LaO .5F .5BiS 2 . Figure H shows x"M at tnree dif- 
ferent temperatures spanning T c = 10.8 K. Our experi- 
mental data show that x (^) does not change within the 
experimental errors when the system transits from the 
normal to SC state. 

In summary, our crystal structure refinements show 
that upon F-doping the buckling of the BiS2 plane re- 
mains very similar to the parent compound. This is in- 
consistent with a theoretical study that the buckling de- 
creases upon F-doping, the change of which facilitates an 
electronic environment for the appearance of the super- 
conductivity [16]. Furthermore, it was theoretically pre- 
dicted that, in the SC phase, a significant change in the 
phonon density of states at low energies would occur due 
to a possible large e-ph coupling. Our inelastic neutron 
scattering data, however, yield no considerable change in 
the low energy phonon modes as the system becomes SC 
either by F-doping or by cooling through the supercon- 
ducting transition. Our results should provide important 
constraints on future theoretical works examining these 
new Bi-based superconductors. 
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